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Abstract 

Background: The invasion-metastasis cascade of cancer involves a process of parallel 
progression. A biological interface (module) in which cells is linked with ECM 
(extracellular matrix) by CAMs (cell adhesion molecules) has been proposed as a tool 
for tracing cancer spatiotemporal dynamics. 

Methods: A mathematical model was established to simulate cancer cell migration. 
Human uterine leiomyoma specimens, in vitro cell migration assay, quantitative real- 
time PCR, western blotting, dynamic viscosity, and an in vivo C57BL6 mouse model 
were used to verify the predictive findings of our model. 

Results: The return to origin probability (RTOP) and its related CAM expression ratio 
in tumors, so-called "tumor self-seeding", gradually decreased with increased tumor 
size, and approached the 3D Polya random walk constant (0.340537) in a periodic 
structure. The biphasic pattern of cancer cell migration revealed that cancer cells 
initially grew together and subsequently began spreading. A higher viscosity of fillers 
applied to the cancer surface was associated with a significantly greater inhibitory 
effect on cancer migration, in accordance with the Stokes-Einstein equation. 

Conclusion: The positional probability and cell-CAM-ECM interface (module) in the 
fractal framework helped us decipher cancer spatiotemporal dynamics; in addition 
we modeled the methods of cancer control by manipulating the microenvironment 
plasticity or inhibiting the CAM expression to the Polya random walk, Polya constant. 
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Background 

Cancer behavior is malignant when the cancer invades neighboring tissues or migrates 
to survive in the ectopic site (metastasis) [1]. The main kinetic features that distinguish 
malignant cancer from benign tumors or normal tissues are that cancer cells do not re- 
main in situ in a stable condition, but spread rapidly into neighboring and distant sites 
through lymphatic or vascular routes. The growth, invasion, and metastases occur con- 
currently in the microenvironment. An emerging paradigm shift challenges the trad- 
itional invasion-metastasis cascade model [2,3] and proposes instead a "parallel 
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progression" cancer model. In this new concept, cancer cells are thought to move by a 
continuous process of "diffusion" or "random walk", although their natural retraction- 
protrusion crawling motion is slower than molecular diffusion. The tracking of single 
molecules has supported the notion that the sub-diffusion of biological matter in non- 
homogeneous and anisotropic microenvironments is non-Fickian [4], 

The context-dependent responses of biological organisms to their microenviron- 
ments, including the ECM (extracellular matrix), generate a continuous sequence of 
states. The dynamic interplay among internal signaling pathways, surface topology, and 
external microenvironments renders biological processes nonlinear. EMT (epithelial to 
mesenchymal transition) and its reverse process MET (mesenchymal to epithelial tran- 
sition) are known to play a role in cancer progression [5], The biological cell-CAM (cell 
adhesion molecule)-ECM interface affects function, and vice versa. In other words, 
each functional cell is related to a set of corresponding molecular complexes in the spe- 
cific biological module [6]. For example, CAMs provoke homotypic connections among 
cells (such as E-cadherins) or heterotypic linkage between cells and ECM (such as in- 
tegrin a v ) in the chemotactic or haptotactic process of cancer cells [7]. Integrins are 
heterodimeric, transmembrane glycoproteins comprising an a chain and a fi chain; at 
least 18 a and 8 /? chains have been identified to date. Integrins represent a key element 
of the biological landscape or surface topology. Integrins also mediate bidirectional sig- 
naling, including the promotion of cell proliferation, adhesion, and resistance to apop- 
tosis. Our previous studies verified that the surface topology of various tissues under 
surface electron microscopy (SEM) was correlated with intracellular complex signaling 
and CAM expressions (e.g., integrin a v ) [8,9]. The concept of the "epigenetic landscape" 
identified by Waddington entails describing perturbed intrinsic phenomena from the 
perspective of the outermost biological surface [10]. Thus, differences in cell-CAM- 
ECM modules wired together correspond to different biological states of cancer. We 
investigated how cancer cells migrate and interact with ECMs in the cell-CAM-ECM 
interface to generate cancer spatiotemporal dynamics [11]. 

Biomathematical models use discrete, continuous, or hybrid approaches to depict 
cancer cell migration behavior, especially the diffusion or random walk mode [12,13]. 
Our model was based on a modified anomalous diffusion theory of cancer cell move- 
ment. In summary, our mathematical model included 2 main points, as follows: 

Fractal space-time scaling of diffusion in probability theory 

The use of probability theory by Albert Einstein (1905) to determine Brownian motion 
provided the first calculation of the Avagadro number [14]. We introduced probability 
theory to provide probability data on cancer cells in different locations at different 
times. There are theorems of probability such as the return to origin probability 
(RTOP), which can be applied to analyze the dynamics of specific cancer cell move- 
ments. Moreover, the probabilistic data of cancer cell movements can be normalized to 
the molecular expressions. Thus, the cell-CAM-ECM interface and its related CAM 
expressions reflect information on time-series alterations of cancer cell position. 

Simple diffusion cannot describe cancer dynamics because real world rarely meets 
ideal conditions (such as symmetry, isotropic or periodic paths of movement, 
homogenous microenvironment resistance, or consistent internal movement drive). For 
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example, the clinical diagnosis of skin cancers has for a long time relied on the "AB" 
acronym, where A refers to asymmetry and B refers to the border (irregular). Uneven 
cell growth and migration are collectively referred to as morphological instability; this 
condition is the hallmark of cancer, both clinically and theoretically [15]. A Darwinian 
clonal selection theory of cancers has verified that the cancerous state represents 
greater heterogeneity and instability than the normal state [16]. Heterogeneity occurs 
along various biological or spatiotemporal scales, and induces spatial variance in cancer 
morphology [17,18]. 

The assumption of even space-time intervals in simple diffusion is inappropriate. Bio- 
logical phenomena occur in space-time, the dimensions of which cannot be represented 
solely by the integers of Euclidean space. Non-linearity can be transformed to a scaling 
problem by space-time fractal factors and diffusion coefficients of specific directionality. 
That is, cancer cells migrate by random walk with various degrees of space-time jump or 
anomalous diffusion [19]. The fractal space-time scaling factor is introduced to quantify 
the extent of heterogeneity or anisotropy. 

The cell-CAM-ECM interface (module) with a key measurable variable 

Because the polarity and diversity of cancer morphological progression cannot be 
traced easily, our strategy is to monitor one of CAMs as a representative factor of can- 
cer that interacts with the specific elements within the microenvironment through the 
biological cell-CAM-ECM interface. The expressions of integrins, which are located in 
areas of surface irregularities, are altered as the surface of the cancer cell changes. Be- 
cause integrins connect cells at distinctive sites in the ECM, different isoforms of integ- 
rins can be used to monitor the changes in the specific cell-CAM-ECM interface that 
correspond to different points in the progression of cancer. For example, integrin a v 
represents a unique feature in the cell-CAM-ECM interface through binding of fibro- 
nectin or vitronectin, which appears as cancer cells attach to the ECM during the inva- 
sion of vessels. This concept originates from the module. We selected key measurable 
variables (such as positional probability and CAM expression) to simplify the analysis. 

Previous cancer research has focused largely on gene regulation and signaling, 
rarely tracing the migration routes or spatiotemporal dynamics in the microenvir- 
onment cues. We built a mathematical anomalous diffusion model in the modular 
and fractal scaling manners. The experiments were also performed to verify our 
model's in-silico predictions orderly. The essential data were simplified to a single 
equation with minimal measurable parameters. Furthermore, the microenviron- 
ment plasticity was manipulated to control the cell-CAM-ECM interface (mod- 
ule). Our findings may shed light on novel cancer treatments. 

Main text 

The analysis of complex dynamic phenomena such as cancer cell migration through 
the cell-CAM-ECM interface (module) requires a simple but efficient methodology. 
Phenomenological models have fewer variables compared with molecular descriptions, 
enhancing their applicability for experimental data. 

Cancer cells migrate in anisotropic and heterogeneous microenvironments. 
Thus, a simple 3D fractal anomalous diffusion equation is built to describe the 



Tsai et al. Theoretical Biology and Medical Modelling 2012, 9:36 
http://www.tbiomed.eom/content/9/1/36 



Page 4 of 12 



dynamics using the Hausdorff derivative. The equation at the cell-CAM-ECM 
interface r{t) gives rise to 

dp ^^ =DA P (r,i) onr(t), 

where p(~f, t) is the positional probability of cancer cells moving in the fractal 
space-time frame (~r = (xi^ l ,X2^ 2 ,x^ 3 ),i = t"), ji and a are fractal space and time 

/D u 0 0 \ 

parameters respectively; D = 0 D22 0 is a diagonal 3x3 diffusion 

V 0 0 D 33 / 

matrix (with D u representing diagonal entries), and A is the Laplacian operator 
[20]. The probability density function of cancer provides data on cancer cell loca- 
tions in the time series. 

The fractal time-scaling transform was proposed by Hoffmann and Ord and is re- 
ferred to as the "internal clock" [21,22]. Assuming the Cauchy problem (Dirac delta 
function distribution p(~f,Q) = 5{xi)8{x2)S{x^)), the following solution is obtained [23]: 

1 td- 1 t t 

P{r,t)=————— i e a— (1) 
^det(D)(477t) 2 

Probability is transformed to the molecular concentration by the normalization 
process, as follows: 

P (T; t )=*{?; t )ij4{?: t )*T= S (T; t )=^>), 

where d(T*, ~t) is the number of cancer cells in a specific location and time, 8(T*,i) is 
the cancer cell concentration on 3D grids, and c(T*,ij is the molecular expression in 
which the molecular concentration ratio equals the molecular concentration in the can- 
cer sample compared with the molecular concentration in the normal tissue. The vari- 
able A is a normalization factor that is based on the number of CAM subtypes that 
occupy each cancer cell surface. Thus, the probability and related CAM expressions, 
such as integrins, can be described in the same manner. 

In summary, this model makes the following novel predictions, which were verified 
by our in silico simulation, in vitro cell, and in vivo mice experiments: 

1. Based on probability theory, we introduced the key concept of RTOP, as follows: 

RTOP « p(0, t) = 3 (2) 

^/det{D) (4ni)~ 2 



Cancer invasiveness does not cause the patient's death, but metastasis to the vital organs 
does. Unsteady tumor cells, whether benign or malignant, may move in several ways 
(growth, invasion, or metastasis). RTOP denotes that the benign or malignant cells migrate 
but some either return to the original location, or at least remain in the vicinity of the ori- 
ginal location. This demonstrated so-called "tumor self-seeding" phenomenon. We focused 
on the garrisoned cells, those that remain in the original location. By contrast, healthy host 
cells remain at the original site and never migrate. RTOP predicts that the CAM 
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concentration ratio c(7^ i) in tumors where a number of migrating cells stay in the vicinity 
gradually approaches a low level as the tumor grows or as time passes. This is depicted by 

c ~ »(0. t a ) = 1 1 . The CAM expressions of tumor surface in the different 

V ' ' y/ AA(D)(*"t?)1 ) 

sampling locations are different because of the different return-to-origin paths and diffusion 
coefficients. However, the data may be constant if the microenvironment with related diffu- 
sion is isotropic and homogeneous. Our finding contradicts the traditional view and the 
tumor-node-metastasis (TNM) cancer staging system, in which molecular expressions in- 
crease as a tumor grows. 

Polya and other researchers have proven that RTOP on a 3D periodic lattice is 0.340537 
[24]. Because the muscle has a relatively periodic structure, we use benign smooth muscle 
tumors to verify whether their CAM expressions are close to the Polya random walk con- 
stant. Specimens of human uterine leiomyoma surface (n = 10) and adjacent normal tissues 
were obtained from patients undergoing resection at Taipei Medical University Hospital 
(TMUH). The study protocol was approved by the institutional review board of TMUH 
(OT-03-08-06). We conducted mRNA and western blot protein analysis to quantify integrin 
av expressions. Scanning electron microscopy (SEM) was further used to examine the sur- 
face of benign tumors. 

2. The transition of spatiotemporal probability provides evidence of alterations in the 
cancer state. Thus, the curve of probability or related cell-CAM-ECM modular 
expressions can show the spatiotemporal pattern of cancer cell migration. 

We cultured three different cell lines (Lewis lung cancer (LLC), human squamous cell 
carcinoma (SAS) and human pancreatic carcinoma- 1 (PANC-1) in the Dulbecco modi- 
fied Eagle medium (DMEM) supplemented with 7.5% fetal bovine serum (FBS). In vitro 
cell migration assay was performed using a BD Falcon cell culture insert (BD Bios- 
ciences). We seeded 1 x 10 5 cells into the upper part of each chamber. Cells on the re- 
verse side of the membrane were stained with 0.1% crystal violet. The migrating cells 
were counted under a microscope at 100 X at 6, 18, 24, 32 and 36 h. 

We also conducted in vivo experiments, for which C57BL/6 mice were purchased 
from BioLASCO Co., Ltd. (Taiwan). When the mice were 6 to 8 wk of age, they were 
injected subcutaneously with 1 x 10 5 LLC cells in the right flank. Our animal study was 
approved under the guidelines established by the Taipei Medical University Ethical 
Committee for Laboratory Animals. Tumor weights were documented on Days 15, 25, 
40, and 47, or on the day of death. 

3. The diffusion coefficient (D) was inversely related to the viscosity {rf) of the 
microenvironment, in accordance with the Stokes-Einstein equation D = (where T is 
absolute temperature; 7 is viscosity; r is the radius of the particle; and k B is the Boltzmann 
constant). Thus, the extent of cancer cell movement depends on the viscosity of 
surrounding microenvironments. Fillers form a barrier to hold down cancer cells. Cancer 
cells must overcome these obstacles by creating the "path-generating" way. 

For the in vitro and in vivo experiments on manipulation of the microenvironment 
properties, mice were randomly divided into 3 groups (n = 5 animals/group for each 



Tsai et al. Theoretical Biology and Medical Modelling 2012, 9:36 
http://www.tbiomed.eom/content/9/1/36 



Page 6 of 12 



time point). The 3 groups received the following treatments: cells only; cells plus HA 
(Restylane Perlane®, Q-Med AB, Sweden); and cells plus polyacrylamide (PAM) 
(Aquamid®, CONTURA, Denmark). The dynamic viscosity of the LLC, LLC + HA 
and LLC + PAM cancer samples (0.5 gm) at different time points was measured using 
a Cone & Plate CAP-2000 Viscometer (Brookfield Engineering Labs., Inc., Middle- 
boro, MA, USA) at 37±0.1°C. The shear rate was 237.6 s" . The unit of dynamic vis- 
cosity is the poise (P). Different fillers were injected around the tumor surface on the 
seventh day of injections. 

Results 

CAM expressions of benign muscle tumors decrease as RTOP simulation 

The simulation of Eq. (2) showed that RTOP (CAM concentration ratio cCr^i)) grad- 
ually approached a low level as the tumor grew or with the passage of time under a 
tumor-specific fractal speed (Figure 1). 

The integrin a v of uterine leiomyoma (benign muscle tumor) was normalized to 
GAPDH and healthy tissues. The values for this parameter showed a descending slope 
and approached a low level as the tumor size increased (e.g., uterine leiomyoma: 
mRNA = 0.39 ± 0.06 and protein = 0.47 ± 0.08) (Figure 2). 

The positional probability (CAM expression) of different cancer cells showed a biphasic 
pattern 

The simulated curve of positional probability and molecular expressions in Eq. (1) 
revealed a biphasic pattern of an initial increase followed by later decline, also known 
as Levy flight (Figure 3). In vitro cell migration assay revealed that the number of three 
different stained cell lines lodged at specific locations underneath the membranes dis- 
played a biphasic pattern, regardless of the type of filler added to the microenvironment 




200 400 600 800 1000 

t(day) 

Figure 1 The simulation indicated that RTOP (return to origin probability) gradually approached as 
low level as the tumor increased (x, =0.5 to 0.7 mm) or as time passed (0 to 1000 days) (D= 10~ 6 / 
d mm 2 ). The fractal factors were a = 0.4 and /3 1|2 ,3 = 0.9. 
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Figure 2 Integrin a v mRNA and protein of uterine leiomyoma (n = 10), which were normalized to 
GAPDH and normal uterine tissues, showed a descending slope and approached a low constant 
level as the tumor size increased (mRNA = 0.39 ±0.06 and protein = 0.47 ±0.08). The experiments 
were performed in duplicate. The surface topology of uterine leiomyomas revealed a relatively smoother 
contour correlated with their gross morphology than normal tissues (inset pictures). 



(Figure 4). In vivo experiment with C57BL/6 mice subcutaneously injected with LLC, 
integrin a v expressions also confirmed this trend (Figure 5). 

Manipulating microenvironment plasticity altered cancer behavior 

The experiments on microenvironment plasticity, together with Eq. (1), showed that 
cancer cell migration depended on heterogeneity or on the diffusion coefficients of 
microenvironments. The Stokes-Einstein equation indicates that higher viscosity (het- 
erogeneity) of materials reduces the diffusion coefficient (p < 0.01, Figure 6). In vivo 
experiments, the fillers caused tumor weight to decrease markedly (p < 0.01, Figure 7). 




Figure 3 In the simulation, positional probability showed a biphasic pattern with an initial increase 
followed by a decrease, for time and location (x, =x 2 = x 3l D r7 = D 2 2 = D S3 = 0.01/dmm 2 , a = 0.4 and 
/3i,2,3 = 0.9). The solid line indicates the curve under a specific Xj =x 2 = x 3 = 05 location value. 
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Figure 4 The numbers of LLC (Lewis Lung Cancer), human squamous cell carcinoma (SAS), and 
human pancreatic carcinoma-1 (PANC-1) cells penetrating the membrane at specific locations 
represented the biphasic pattern in the cell migration assay. The pictures show LLC, SAS and PANC- 
cells at 32 hours in each corresponding figure. All data represent the means ± standard error (S.E.) from 
duplicate experiments (n = 3 for each time point). 



Statistical analyses 

Statistical results are presented as the mean ± standard error (S.E.). To compare the 
data from individual cancer groups, differences in tumor weight were tested by the 
Welch test (overall comparison) and post hoc Dunnett test (multiple comparisons). Dif- 
ferences in dynamic viscosity among the filler groups were assessed using the non- 
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parametric Kruskal-Wallis H test and post hoc Dunnett test. All tests were 2-tailed and 
differences were considered significant if p < 0.01. 



Discussion 

Tumorigenesis involves multiple and complex network mechanisms. To date, single 
targeted molecules have not yielded the expected rewards for diagnosis and treatment. 
Our modeling strategy simulated the so-called drunken unstable cancer cell migration 
as a fractal diffusion or random walk. Fractal factors can quantify the heterogeneity and 
anisotropy of tumors and their microenvironments. The positional probability can eas- 
ily be normalized to a surface molecule expression (e.g., CAMs) of cancer cells as a 
cell-CAM-ECM module because cells and microenvironments reshape and interact 
with each other. Three possibilities describe the location of cancer cells, namely (1) stay 
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in situ and coalesce to a mass; (2) move but return to origin; and (3) migrate to neigh- 
boring or distant sites. Briefly, we posit a new cancer anomalous diffusion hypothesis, 
as follows: motility dynamics underlie a cell-CAM-ECM module and relevant malig- 
nant/benign behaviors. 

Researchers have long understood that dynamical reciprocity (DR) creates various dy- 
namic space-time patterns. The monitoring of CAM expressions seems to provide net 
biological information without myriad complex intracellular signaling detections or tar- 
geted treatments. If no cancer cell exists or lodges in a specific space-time, no related 
CAM is generated to adhere to ECM. This phenomenon is accepted in anoikis (integ- 
rin-mediated death and caspase-dependent apoptosis) [25]. Namely, a specific CAM 
level corresponds to a unique cell-CAM-ECM interface (module). 

Our model simplifies the construction and obtains profound results based on these 
assumptions. First, the concept of RTOP has been well applied to explain why the fila- 
ments in the molecular motor ratchet can always return to their original position after 
contraction or action in one or two dimensions [26]. In our study, RTOP is applied to 
explain the probability of benign or malignant cells which stay, seed or return to the 
original position or solid mass state (tumor self-seeding). Benign tumors always hold 
their solid tumor shape; but the shape or margin of malignant tumors (cancers) blur. 
Several benign tumors have self-limiting invasion features, including uterine leio- 
myoma, breast fibroadenoma, benign prostate hypertrophy, and brain meningioma. 
Namely, CAM expressions in the benign tumors steadily follow our RTOP prediction 
compared to the malignant tumors. The simulation revealed that CAM concentration 
ratio decreased as a benign tumor enlarged in the random structures of the microenvir- 
onments. If we selected a specific microenvironments (the periodic lattice) from all 
possible biological microenvironments, the data of benign tumors (uterine leiomyoma) 
growing in the periodic muscle structures approached the 3D Polya random walk con- 
stant (0.340537). The RTOP value may be clinically significant to the treatment goals 
of molecular target inhibition. Namely, the tumor state may be suppressed to a state as 
characterized by the CAM expressions with little spatiotemporal variance in a regular- 
structure microenvironment. 
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Second, the biphasic pattern of the positional probability indicates the main process 
of cancer progression, namely initial growth with later spread, with a continued parallel 
cancer progression. Cancer cells are initially located in certain vicinity with a peak at a 
specific space-time; subsequently, most cancer cells begin to escape. Non-linearity vio- 
lates chronological sequences of events; however, the main pattern remains evident. 
CAM expressions of cancers may help indicate the state and pattern of cancer in the 
clinical implications. 

Third, the Stokes-Einstein equation provides an evaluation guide to microenviron- 
ment plasticity or the use of applied fillers in cancer control. Artificial control of micro- 
environments, including boundary conditions in differential equations, can impede the 
spread or migration of cancer cells. Our experiments showed that higher viscosity of 
fillers was associated with slower cancer migration or tumor growth. Manipulating 
ECM components provides information on how cancer progression might be attenu- 
ated. The results of our modeling experiments, especially those for RTOP and the use 
of fillers, indicate that our mathematical model may be useful in the development of 
new strategies for cancer prevention and treatment. 

Conclusions 

We propose that complex diseases require dynamic, modular assessment and control. 
Biological networks appear to exhibit modularity as topological structures. The inter- 
play between the positional probability and cell-CAM-ECM module shown in our study 
led us to establish a simplified mathematical model, and further provides useful infor- 
mation on cancer progression. 
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